
THERMAL CONTRACTION 
OF MERCURY 

by Han-Shou Liu 

Goddard Space Flight Center 
Greenbelt, Md, 20771 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • JUNE 1970 


TECH LIBRARY KAFB, NM 


1 . Report No. 


2. Government Accession No. 


NASA TN D-5867 1 

4. Title and Subtitle 

Thermal Contraction of Mercury 


7. Author(s) 

Han-Shou Liu 

9. Performing Organization Name and Address 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 


12. Sponsoring Agency Name and Address 


National Aeronautics and Space Administration 
Washington , D. C . 20546 

15. Supplementary Notes 


16. Abstract 


TECH LIBRARY KAFB, NM 




3. Recipicm s catalog No. 

5. Report Date 

June 1970 

6. Performing Organization Code 

8. Performing Organization Report No. 

10. Work Unit No. 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

14. Sponsoring Agency Code 



The cooling of the planet Mercury would produce a secular decrease of the moment 
of inertia and accelerate its rotation. It is shown that the initiation of the tidal evolu- 
tion on Mercury was balanced with the effect of the thermal contraction during the 
earlier epoch of cooling in a period of about 48 million years. The theory also serves 
to explain the surface features of Mercury. 


|17. Key Words Suggested by Author 


18. Distribution Statement 


Unclassified-Unlimited 


19. Security Classif. (of this report) 

. . j 

20. Security Classif. (of this page) 

21 . No. of Pages 

Unclassified 

l Unclassified 

9 


22. Price* 

$3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 

Springfield, Virginia 22151 




CONTENTS 


Page 


INTRODUCTION 1 

COOLING 1 

SECULAR DECREASE IN MOMENT OF INERTIA 2 

AXIAL ROTATION 4 

OROGENETIC SURFACE FEATURES 5 

CONCLUSIONS 8 

ACKNOWLEDGMENT 8 

References 8 


Appendix A-List of Symbols 


9 



THERMAL CONTRACTION OF MERCURY 


by 

Han-Shou Liu 

Goddard Space Flight Center 


INTRODUCTION 

Liu (1968) has suggested that the trapping of the rotation of the planet Mercury was originally effected 
by thermal expansion or contraction of its configuration during solidification. According to this hypothesis, 
Mercury was orignally at a very high temperature, and different parts have since cooled by different amounts. 
This suggestion was based on results of quantitative calculations. In this paper, the derivation of the re- 
sults on the underlying physical problem is given. On the basis of the analysis made in the course of this 
study, it is shown that the thermal process on the planet Mercury during the earlier epoch of cooling is quite 
adequate to account for the origin of trapping. The theory also predicts the existence of folded and thrusted 
mountain systems on the surface of Mercury. 

COOLING 

The thermal history of Mercury, after solidification, can be determined from the laws of thermal conduc- 
tion. The equation of conduction is 


dT Id 
P ° dt ~ r 2 dr 



4 - A (f, t), 


( 1 ) 


where p is the density, c is the specific heat, T is the temperature, t is the time since solidification, r is 
the radial distance from the center of the planet, k is the thermal conductivity, and A is the rate of genera- 
tion of heat per unit volume. 

The general solution of Equation 1 may be separated into two parts. One part concerns a nonradio- 
active Mercury that cools by conduction from the actual initial state; the other concerns a Mercury initially 
at zero temperature which is warmed by radioactivity and loses heat by conduction to the surface. Numerous 
solutions can be given. For the present purposes it is convenient to make certain approximations which can- 
not introduce much error and which lead to an easily intelligible result. One may assume that (1) p, c, and 
k are the same at all depths; (2) cooling is appreciable only at depths which are a small fraction of the radius 
of Mercury; (3) A is independent of time and is appreciable only in a layer whose thickness is very small com- 
pared with the depth to which cooling has extended. With these simplifications, the temperature below the 
radioactive layer at a depth x is (Jeffreys, 1959) 



(?) 


T = [iX + 



where p is the initial temperature gradient, S is the initial surface temperature, 1 is the depth of radioactive 
layer, h 2 = k/ pc is the thermometric conductivity, t is the time since solidification, and the error function is 
defined by 


rO 

erf 6 =— / exp (-y 2 ) d y. 

7T h JO 


For 


A = 7.3 x 10" 13 cal cm' 3 sec" 1 
c = 0.2 cal g" 1 deg' 1 
k = 0.5 x 10' 2 cal cm' 1 sec" 1 deg" 1 
1 = 1.2 x 10 6 cm 
S = 1400° K 
p = 5 g cm" 3 

it is found that 

T = px + L erf \ + Q 

\2/it /2 / 

and the rate of cooling is 

k. — 2 — exp l-± V 

(9t 77^2 2At3/2 ^ 4/,2( ) 

where 

h 2 = 5 x 10" 3 cm 2 sec' 1 
L = 1295° K 
Q = 105° K 

p = 4 x 10" 5 deg cm" 1 . 


( 3 ) 


(4) 


SECULAR DECREASE IN MOMENT OF INERTIA 

Throughout the process of cooling, a certain amount of mechanical adjustment must have occurred in the 
shape of Mercury. Consider a shell of internal radius r and external radius r + dr. The rise of temperature 
of the shell in a definite time is of course a very small negative value. The density of the shell will change 
at the same time from the initial value p to p (1 - 3a T) where a is the coefficient of linear thermal expansion. 
If the inner radius changes to r(l + /3), the external radius becomes 

r (1 + f$) + dr |l + 

Hence the mass of the shell after the change of temperature is 


2 



( 5 ) 


477f2 p dr |^1 + 2/3 +-J-(r/3) - 3aT| 

But the mass is unaltered. Thus the equation of continuity is 

2/3 + j r (r/S) - 3a T = 0. 

For a given cooling, this is the equation for determining the value of /3 subject to the boundary condition 
that j3 is zero at the center of Mercury. Therefore, for any shell, (3 is determined by the changes of tempera- 
ture within that shell. If the shell simply expands independently of the interior, the radius will increase by 
aTr instead of by r/3, so that the excess 


r (/3 - aT) = dr 

is due to stretching. The stretching is given by (Jeffreys, 1959) 

-- - (8r3) = r^^H. 
dr dr 

The rate of stretching below the radioactive layer at a depth x is then 



where r Q is the radius of Mercury. Now, the value a can be expressed as 

° - e i + f 2 ^ 

where € l and e 2 are constants. From Equations 3, 4, 5, and 6, one obtains the amount of stretching at the 
surface 

8 = - ^ | ! 0rt) I/2 + f 2 [(Q + 2 1/2 L) Grt)* + ixhnt^. (7) 


For 


fl = 7x 10' 6 deg' 1 , f 2 = 2.4 x 10' 8 deg' 2 , and t = 1.5 x 10 15 sec, 


Equation 7 gives 


5 = - 2.8 x 1CT 3 . • 


From this the increase of the radius of Mercury by thermal expansion during the earlier stage of cooling of 
the shell with a depth of 300 km, is at once found to be about -7 km. The surface of Mercury, by crumpling, 
has therefore been diminished by about 0.2 x 10 6 km 2 . The rate of increase of the radius is, then, about 
— 4.7 x 10“ 10 cm sec' 1 . Therefore the average rate of change of the moment of inertia about the axis of rota- 
tion during the cooling period of about 48 million years is approximately 


3 



= - 3.1 x 10" 10 M m r Q cm 2 g sec" 1 , 


( 8 ) 


where M m is the mass of Mercury. 


AXIAL ROTATION 

During the earlier epoch of cooling, Mercury would rotate about its polar axis with a time-dependent in- 
ertial tensor. The rotation of Mercury is described by 


d Tr d(/ + <w l m m 

dT c <«> — dt — = ~ N (9) 

— 

where / is the true anomaly, 4> is the angle of orientation of Mercury relative to the Sun, and N is the couple 
due to solar tides. Equation 9 can be rewritten in the following form 


d 2 (/ + 9 4) N _H_ 

dt 2 C (t) C (t) 


( 10 ) 


where 


18 M/ 

N=^-k 2 nG p — 
52 M ‘ 


r 0 6 (1 + e cos f) e 
a 6 (1 - e 2 ) 6 


‘ C (tV sin 2<r 


and 




d (/ + *) dC 


it) 


dt 


dt 


In Equation 10, a is the major semiaxis, e is the orbital eccentricity, G is the gravitational constant, k 2 is 
the Love number, € is the phase lag of the tides, and M s is the mass of the Sun. 


For e = 0.2; G = 6.7 x 10 8 dyne cm 2 g' 2 ; k 2 = 0.02; 2e = 0.005; p = 5.0 g cm" 3 ; M s /M m = 6.0 x 10 6 ; 
r Q /a =4.1x 10" 5 ; and r Q = 2.42 x 10 8 cm, the average value of W/C (jt) over a period of revolution is 

=0 (10“ 23 ) rad sec' 2 . 

The calculation also includes the reversal of the phase lag during perihelion passage* (Peale and Gold, 
1965; Liu, 1966). 

At the 3:2 resonance condition, d(/ + c/>)/dt = (3/2)n = 1.8 x 10' 6 rad sec" 1 , where n is the mean orbital 
motion. From Equation 8 one obtains 


*Liu, H-S., unpublished results, 1966. 
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= - 0 (10 23 ) rad sec' 2 . 

C U) 

Therefore, it is seen that the initiation of the tidal evolution was completely balanced by the effect of the 
thermal contraction during the epoch of cooling lasting about 48 million years. 

In applying Darwin's hypothesis of the formation of the Earth, it should be noted that in the initial stage 
the thermal contraction of Mercury may be more effective in determining its rotation than the solar tides be- 
cause the loss of heat at that stage was exceptionally great. Therefore the estimate of H/C^ may be re- 
garded as a lower limit. Since the calculation for this quantity is of the same order of magnitude as that for 
N/C( t ), it seems reasonable to conclude that the influence of thermal contraction provided the conditions 
for the rotation of Mercury to be trapped into the 3:2 resonance state. 

OROGENETIC SURFACE FEATURES 

By performing the integrations on the right-hand side of Equation 6, one finds 

^=U l + t 2 Q + 2<r 2 fiht /2 q + e 2 L erf q) exp (- q 2 ) 

- ) Oj + ( 2 Q) exp (— < 2 2 ) 

0rt)*r o \ 1 2 

+ (j2qf exp (- q 2 ) + n k (1- erf q)J 

+ e 2 L^exp (-q 2 ) erf q + 2 V2 (1-erf 2 V2 q)jj,, (11) 

where 


q = 


x 

2ht /2 ' 


The value of 38/ dt is negative when q is zero, but when q is large enough, 38/ 3t is positive. For one value 
of g, 38/ dt is zero, and the material is being neither crumpled nor stretched. The level with this value of q 
is defined as the level of no strain. For 38/3t = 0, Equation 11 gives q = 0.1. Hence, the depth of the level 
of no strain is about 40 km below the surface at the present time. 

The internal states of strain, which were developed by thermal contraction, are illustrated in Figure 1. 
The mass from the center of Mercury to within about 300 km below its surface has not been affected by any 
appreciable cooling. Within the shell from 300 to 40 km below the surface, cooling by conduction is taking 
place, and this shell is contracting about the unchanging interior. Hence, it is in a state of internal tension. 
Above the level of no strain, the layers have already cooled and therefore are under a state of internal com- 
pression with a horizontal crushing stress. Actually, in a solid material with a finite strength, differences 
of cooling would produce stress-stresses which would increase until they reached the strength of the mate- 
rial and then be mainly relieved by fracture. The power freed by the fracture of the contracting solid shell 
is considered here to be the cause of Mercurian orogenesis. 
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Figure 1 - States of thermal strain of Mercury. 

Since the solid contracting shell has to fit the unchanging interior, it is necessary to determine whether 
a readjustment that will reduce all great circles of the sphere by different amounts is possible. This can be 
done. For simplicity, suppose that there is crumpling across the meridian of the thermal bulge and also 
across the equator from 90° E through 180° E to 90° W. Every great circle of the sphere would intersect both 
of these arcs. Because of the solar heating after the rotation becomes commensurable with the mean motion 
(Liu, 1968), the shortenings of different great circles would vary at different positions. Such a distribution 
is remarkably like that shown in Schiaparelli’s map of Mercury (Figure 2). (Schiaparelli, 1889). 

What is likely to happen when a fracture is formed at the surface by contraction? The stress differences 
arise because the outer shell is too large to fit the interior; they can be relieved only by a displacement that 
reduces the outer surface. A vertical fracture or a horizontal split will not do. There must be an oblique 
fracture with sliding along the thrust plane. Relief by thrusting would be at an angle of somewhat less than 
45° to the horizontal because the outer shell is in a state of compression (Figure 3). In this way the crust of 
Mercury to a depth of about 40 km would be overthrusted, leaving the upper edge overhanging by a similar 
height. Therefore the surface features of Mercury (Figure 2) are probably folded and thrusted mountain sys- 
tems that were formed as a natural consequence of thermal contraction. 
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CONCLUSIONS 


The version of the thermal contraction hypothesis presented in this paper appears to offer a sound 
physical theory with which predictions for Mercury can be made and checked. 

ACKNOWLEDGMENT 

The author expresses his appreciation to Dr. John A. O'Keefe for his comments. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, December 18, 1969 
188-48-01-05-51 


REFERENCES 

Jeffreys, H. , The Earth (New York: Cambridge University Press, 1959). 

Liu, H-S., “Mercury Has Two Permanent Thermal Bulges/' Science, 159: 306, 1968. 

Peale, S. J., and Gold, T., “Rotation of the Planet Mercury," Nature , 206: 1240, 1965. 
Schiaparelli, G. V., “Sulla Rotazione Di Mercurio," Astronomishe Nachrichten, 123: 2944, 1889. 


8 



Appendix A 


LIST OF SYMBOLS 


= major semiaxis 

= rate of generation of heat per unit volume 

= specific heat of body 

- moment of inertia about axis of rotation 

= orbital eccentricity 

= true anomaly 

= gravitational constant 

= thermometric conductivity, — 

pc 

= thermal conductivity 
= Love number 
= depth of radioactive layer 

= temperature difference between surfaces of radioactive layers 
= mass of Mercury 
= mass of Sun 
= mean orbital motion 
= couple due to solar tides 
= temperature change across radioactive layer 

= internal radial distance of a shell from center of plane at which temperature is determined 
= radius of planet Mercury 
= initial surface temperature 
= time since solidification 
= temperature of body 
= depth below surface of radioactive layer 
= variable of error function 
= coefficient of linear thermal expansion 

= amount of stretching as a function of time and the depth of the planet 
= phase lag of the tides 
= constant (value of a at T = 0) 

= constant (change with T in value of a) 

= limit of error function 
= initial temperature gradient 
= density of mass 

= angle of orientation of Mercury relative to the Sun 
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